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^^ Abstract 

^^ After a brief description of the basic theoretical concepts of the string 

model approach, we discuss the limitations of the method, at low energies, at 
*K^ high energies, and at high particle densities. We also report on recent efforts 

CO to overcome these limitations. 

o 

--H 1 Introduction 

fn During the past decade the string model has been quite successfully applied to 

describe hadronic interactions, in particular nucleus-nucleus collisions at CERN- 
SPS energies. Awaiting experiments with larger nuclei at higher energies in the 
future, it appears useful to reflect about the limitations of the appraoch and about 
the efforts to go beyond. 

In chapter 2, we present some overview of modelling hadronic interactions, also 
to clarify what we actually mean by "string model approach" and to show relations 
to other models. In chapter 3 we discuss the basic concepts of the string model, 
with emphasis on the way the string is introduced, namely as a parametrization of 
squared amplitudes and not as an elementary object. We then discuss in chapter 4 
the limitations of the approach, and recent efforts to go beyond. 
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2 Overview 

When modelling ultrarelativistic liadronic interactions, one first needs to specify 
whether soft or hard interactions are to be considered. "Hard" refers to scatterings 
involving a large momentum transfer - large enough to allow for a perturbative treat- 
ment, otherwise, for small momentum transfer, the term "soft" is used. Whether a 
particle originates from soft or hard scattering, can be most easily seen by inspect- 
ing transverse momentum (pt) spectra. Pion spectra typically show an exponential 
low-pt behaviour (for pt < 2 GeV/c) and a laige-pt behaviour of the form p^"' (see 
fig. 1). Such a power law behaviour is typical for hard processes, described in terms 
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Figure 1: Tranverse momentum distri- 
bution separating soft and hard scatter- 
ing. 



Figure 2: Realative weight of hard 
scattering versus energy. 



of elementary QCD diagrams, whereas the exponential low-pt part is due to soft 
processes. The relative weight of hard processes is energy dependent (see fig. 2): 
at SPS energies (20 GeV) hard scattering plays only a marginal role, at LHC (6300 
GeV) hard scattering is crucial. At RHIC energy (200 GeV), hard scattering is 
important, however, the interaction is still most likely soft. For the following, we 
restrict ourselves to soft interactions, keeping in mind that this provides a descrip- 
tion up to moderately high energies only (< 100 GeV), but needs extension to go 
beyond. 
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Table 1: Classification of models for ultrarelativistic (UR) hadronic interactions 



There are few different approaches to describe soft interactions, hadron-hadron as 
well as hadron-nucleus and nucleus-nucleus scattering (see table 1). Models based on 
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molecular dynamics (MD) are a straightforward generalization from lower energies 
[1, 2]. On the other hand, there are string models, either properly defined in the 
frame work of relativistic quantum theory - based on Gribiv-Regge-Theory (GRT) 
[3, 4, 5, 6, 7] - or string models introduced in a more phenomenological fashion [8]. 
We do not want to comment on MD models nor on phenomenological string models, 
we rather provide a critical discussion of GRT-based string models, its theoretical 
justification and apparent limitations. So in the following, "string model" refers 
solely to GRT-based string models. 



3 The String Model of Hadronic Interactions 

In the following we discuss the theoretical concepts of the string model approach, 
implemented in models like VENUS, DPM, and QGSM. One first considers the 
elastic amplitude A(s,t) for hadron-hadron scattering, with s and t being the Man- 
delstam variables s = [pi +^2)^ and t = [pi — ps)^, where p\,P2 are the incoming and 
P3tP4 the outgoing momenta. The amplitude A(s,t) is given as a multiple Pomeron 
exchange series. 




with the zigzag lines symbolizing Pomerons. The Pomeron is the elementary ex- 
change object here, without being elementary in terms of quarks and gluons. Ori- 
ginally the Pomeron was thought to be a gluon ladder. According to Veneziano, a 
Pomeron is the sum of all QCD diagrams of cylindrical topology A gluon ladder is, 
by the way, of cylindrical topology, in some sense the simplest cylinder. 

Whatever the precise nature of the Pomeron may be: first of all one simply 
parametrizes the Pomeron propagator as 

G'(s,t)~s"W = s"(°)+"'*, (1) 

with this Regge pole form coming from general considerations of amplitudes in 
the limit s — )• 00. Two parameters characterize the Regge trajectory a[t): the 
"intercept" a(0) and the "slope" a', both being adjusted to fit data. With this 
simple form for G, the convolutions G . . . G can be worked out and summed over, 
and one obtains 
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with 

7{s,b) = l-e-(^''\ (3) 



String models ... limitations and new directions 





19.4 GeV 


200 GeV 


6.3 TeV 


Wi 


0.54 


0.52 


0.47 


W2 


0.24 


0.23 


0.22 


W3 


0.12 


0.13 


0.13 


W4 


0.059 


0.067 


0.081 


W5 


0.026 


0.033 


0.047 


We 


0.011 


0.015 


0.026 


W'j 


0.0040 


0.0061 


0.0130 


Wg 


0.0013 


0.0023 


0.0060 


Wg 


0.00042 


0.00080 


0.00260 


Wio 


0.00012 


0.00025 


0.00100 



Table 2: Weights Wm for cutting m Pomerons 
(m colour exchanges). 



where uj is the Fourier transform of the Pomeron propagator G, 



u(s,b) = — I d''kG(s,t)e 

ITT 



-ikh 



(4) 



It should be noted that only two-dimensional integrations occur {(Ph^cPk). The 
reason is that at high energies transferred momenta are purely transverse, and there- 
fore longitudinal and transverse degrees of freedom decouple. Longitudinal integra- 
tions can be performed first, so the convolutions G ® ... ® G mentioned earlier are 
simply two-dimensional integrations referring to transverse momenta. This is the 
reason why results can be presented in "impact parameter representation", namely 
in the form j (Ph.... 

Using the above result for A(s,t), the elastic cross section is given as 

'^"^^ ^d%HsM' (5) 



C^el 



dt- 



dt 



and, using the optical theorem, the total cross section can be written as 

(jtot = 8^ Im A(s,0) = I d^b2Re-f{s,b). 
As a consequence, the inelastic cross section is given as 



C^inel — CTtot — CTel 



d%{l 



-2cu(s,b) 



}• 



(6) 



(7) 



After having discussed elastic scattering, we would like to study inelastic amp- 
litudes 
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representing particle production. Unitarity relates elastic and inelastic amplitudes 
via 



2 Im 



where the sum extends over a complete set of final states. Since inelastic amplitudes 
cannot be calculated directly, one first treats the imaginary part of the elastic amp- 
litude and then draws conclusions about inelastic scattering. For this purpose, 
the AGK technique is employed [9], which provides a method to express ImA(s,t) 
in terms of elementary inelastic scatterings, given as lniG{s^t)^ with G being the 
Pomeron propagator. Using this technique and the optical theorem, the total cross 
section can be expanded as 



Cftot 



E 

?Ti = 



(Tr, 



(8) 



where 0"^ is the cross section for m elementary inelastic scatterings. The 0"^ can be 
calculated, they are weakly energy dependent, as shown in table 2, where we show 



Wr^ 



C^m/cTi 



inel* 



To be more specific about the "elementary inelastic processes", we need to know 
something about the nature of the Pomeron. We adopt Veneziano's picture of a 
Pomeron being a cylinder (of gluons and quark loops). 





As a consequence, the imaginary part of the Pomeron amplitude is taken to be a 
squared "cut cylinder". 



2 Im 



^ 



This is somehow the unitarity equation of a single Pomeron (compare with the full 
unitarity equation). Such a cut cylinder is still too complicated to be evaluated. It is, 
however, well known that planar QCD diagrams have surprisingly similar properties 
as compared to fragmenting strings [3], and correspondingly the cut cylinder is 
parametrized as two independently fragmenting relativistic strings. To indicate this 
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Figure 3: Quark line dia- 
gram. 



and to keep track of the flavour flow, it turned out to be useful to to introduce 
quark line diagrams as shown in flg. 3 for a single colour exchange in the case of the 
VENUS model (Similar diagrams are also used in the DPM and QGSM). Quarks 
(and antiquarks) are represented by horizontal lines, and vertical connection of quark 
lines represent strings (including the incident hadrons). The quarks belonging to a 
string are highlighted by dots. This interaction process represents a colour exchange 
mechanism: before the interaction, the projectile and the target quarks are singlets, 
afterwards - as if colour would be exchanged between a projectile and a target 



Pi 



Ti 



T, 



ffl (> 

O (> 

© • 

© • 

© 7 • 

<B © A 

O O • 

® f t 

© i 

i .. i 

© 4 

o «> 

i ;; i 

© 4 

© i 

© © i 

o «> 

© i 



Figure 4: Quark line diagram for a nuclear interaction. 
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quark - a diquark forms a singlet together with a quark from the other nucleon. 
The colour exchange is indicated by an arrow in the figure. The generalization to 
multiple colour exchange is straightforward and obvious. 

So far we have been discussing hadron-hadron scattering. The generalization 
to nucleus-nucleus is straightforward: starting from an expression for the elastic 
amplitude in terms of multiple Pomeron exchange, one obtains, by using the AGK 
technique, a multiple colour exchange picture for nuclear scattering. No new para- 
meters need to be introduced, apart of those specifying the nuclear geometry. A 
simple quark line diagram for string formation in nuclear collisions is shown in fig. 
4. 



4 Limitations and New Directions 

The model described so far is the so-called "independent string model", referring to 
the fact that the strings - even in nuclear collisions - fragment independent of each 
other. There are a couple of obvious limitations of the approach: 

• The low energy limit: below some energy ^s^i^ the string picture breaks 
down completely. 

• The high energy limit: beyond some energy ^s^ax new components have 
to be added (— )• "soft & hard scattering") 

• The high density limit: for high densities of produced hadrons, the concept 
of independent strings has to be modified (— )• "secondary interactions"). 

In the following, we discuss these limitations as well as new developments to go 
beyond. 

4.1 Breakdown of the string model at low energies 

Let us first introduce some time scales: the "hadronization time" t^ is the time it 
takes before a hadron is formed, the "inter-collision time" Tc is the time between 
successive nucleon-nucleon collisions, and the "traversal time" tr is the time it 
takes for a projectile nucleon to traverse a nucleus of radius R. To be specific and 
to simplify the discussion, we consider in the following a nucleon-nucleus collision, 
in the rest frame of the nucleus. 

We may define an incident energy leading to 



Th < Tc 

as "low energy". This relation implies that, after a nucleon-nucleon interaction, 
hadrons are produced essentially immediately before another target nucleon gets 
involved. This leads to the "cascade picture": the first nucleon-nucleon interaction 
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produces a couple of hadrons, each of which - when hitting another target nucleon - 
will interact and produce hadrons itself, and so on. Such a scenario can be formulated 
completely in terms of hadron, without referring to the hadronic substructure. 
The situation is completely different at high energy, defined via 

Th > Tr. 

This means the following: after the first interaction, the leading "object" - the 
forward end of a string - traverses the whole nucleus, making interactions whenever 
a target nucleon comes in its way, before a leading hadron is formed outside the 
nucleus. This is actually the scenario employed by the string model approach, and 
therefore the relation t^ > tr may be used to estimate the energy Emin below which 
the string model breaks down: 

p fm 

n = ItI = r^* ^ E—— >TR^2R^ fOfm 



which leads to 



or, in the cms. 



E> lOGeV =:E^in , 



2mjv£'inin ~ 5 GeV 

For energies less that that, we also have to deal with very low mass strings, adding 
another uncertainty. A third reason for not to use the string model below s^in is that 
the Gribov-Regge formalism, the basis of the string model, is based upon asymptotic 
considerations for s/GeV >> f . It appears to be most difficult if not impossible to 
fix all these shortcomings, so we consider s^in ~ 5 GeV an absolute lower limit, and 
do not attempt to go beyond. 

4.2 Soft and hard scattering 

From collider experiments at CFRN and Fermilab we know that for energies (v^) 
beyond 10^ GeV a new feature shows up in nucleon-nucleon scattering: jets, refer- 
ring to almost coUinear bunches of hadrons. The origin of this phenomenon are 
elementary interactions of the type parton + parton — )• parton + parton, where 
each outgoing parton decays into a jet of hadrons. The inclusive cross section CTjet 
for jet production can be calculated by evaluating elementary QCD diagrams, and 
the result is - with some uncertainty due to an arbitrary pt-cutoff - the following 
[10]: CTjet is zero up to -^s sa 30 GeV, increasing slowly up to CTjet ~ 5 mb at 10^ 
GeV, and rising quickly with energy, exceeding 40 mb around 10"^ GeV. The total 
cross section, on the other hand, rises slowly, from 40 mb at 30 GeV, via 45 mb at 
10^ GeV, to around 70 mb at 10"^ GeV. The value of CTjet relative to crtot indicates 
the relative strength of the hard processes, leading to jets. For values of s below 
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Figure 5: Hard scattering. Figure 6: Multiple hard scattering. 

30-100 GeV, hard processes can be ignored. This provides an upper energy hmit for 
the string model, representing pure soft interactions: 



^/I^^ 30- 100 GeV 

To treat higher energies, the string model needs to be extended to incorporate 
also hard processes. 

Let us first briefiy explain the standard perturbative QCD (PQCD) treatment 
of hard processes. Starting point is the following expression for the inclusive cross 
section for jet production [10, 11, 12, 13, 14, 15]: 

^i^ = Y. J dxi J dx2 J dt^Mxi,Q')Mx2,Q') (9) 

with the corresponding diagram shown in fig. 5. The indices i and j refer to partons 
of the nucleons (quarks/antiquarks of a certain fiavour or gluons), fi and fj are the 
corresponding momentum distribution functions, with Xi and X2 being the longitud- 
inal momentum fractions with respect to the incident momenta, and with Q being 
the transferred momentum. The quantity daij/dt represents the elementary cross 
section, calculated from summing QCD diagrams of lowest order. 

It turned out that the inclusive cross section at high energies exceeds the inelastic 
cross section, 

C^h > CTinel, (10) 

which refiects multiple scattering: The ratio 

n:=^ (11) 



C^inel 



represents the average number of hard scatterings - which may be larger than one. 
In order to quantitatively reproduce not only pt spectra but also multiplicity dis- 
tributions, an impact parameter model has been introduced [10, 11, 14]. Here one 
introduces an average number of collisions for a given impact parameter as 

n{b) := a^T{b) (12) 

with 



T{b):= p{h-h')p{h% (13) 
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where 

p{b) := / dzp{K,by,z) (14) 



is the transverse matter density of a nucleon. The quantity T[b) is a measure of the 
overlap of the two nucleons: T[b) is largest for 6 = (representing complete overlap), 
drops with increasing 6, and approaches zero around b = 2i?nucieon (representing two 
touching nucleons). 

Assuming independent multiple scattering, one writes the probability prob(j, b) 
of j scatterings at given impact parameter as a Poissonian, 

prob(j,6) = ^'e-'^W, (15) 

with h[b) given in eq. (12) (see fig. 6). Summing over j and integrating over 6, one 
obtains the inelastic cross section, 



/oo 
o?^6^prob(j, 6) 

= I d%{l-exp[-ai,T{b)]} . (16) 

This is the well-known eikonal form. 

We realize a large formal similarity between soft and hard scattering: in both 
cases multiple scattering leads to the "eikonal form" of total or inelastic cross section. 
For soft scattering, we have multiple cut Pomerons (multiple colour exchanges), for 
hard scattering, we have multiple jet production. Combining soft and hard amounts 
to considering a Pomeron (P) having a soft and a hard component, 

P = Psoft + Phard, 

with a hard cut Pomeron representing a hard scattering and parton radiation. There 
are several attempts along theses lines [3, 7, 13], not being satisfactory though due to 
the necessity of an artificial low-pt cut-off. There is at present also a large uncertainty 
concerning predictions for multiplicities in nucleus-nucleus collision at RHIC or LHC 
energies. 

4.3 Secondary Interactions 

It is obvious that the independent string model has to break down, when the density 
of produced hadrons from string decay gets too large. This seems to happen already 
for moderately heavy nuclei [S -\- S) at SPS energies. 

Three approaches have been suggested to extend the independent string model 
correspondingly: 

• the hadronic cascade: 
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• the string fusion mechanism; 

• the formation of quark-matter clusters. 

The hadronic cascade [7, 2] amounts to following the trajectories of hadrons 
from string decay and considering interactions according to tabulated cross sections, 
when two trajectories come sufficiently close. This method of purely hadronic binary 
interactions, however, obviously breaks down for densities such that more than two 
hadrons are likely to be closer than an interaction distance (which is probably the 
case already for S+S at SPS energy). 

String fusion [2, 8, f 6] considers secondary interactions already at an early stage, 
before particles are formed. In case of two strings being close together, they "fuse" 
to form a new string with increased charges at the endpoints. Two 3 — 3 strings 
will fuse into a 6 — 6 string (with weight 2/3) or a 3 — 3 string (with weight 1/3), 
where the 3, 3, 6, and 6 refer to SU(3) representations, and the weights 2/3 and 
1/3 result from the rules of multiplying SU(3) representations (3x3 = 6 + 3). The 
break of a 6 — 6 string requires at least a diquark-antidiquark production (which 
is rare for "normal strings") and will lead to a large baryon-antibaryon production 
rate. Details about string fusion may be found in [2, 8, 16]. 

Formation of quark matter clusters (QM-cluster approach) is the third approach 
to extend the independent string model [17], to be discussed in the following. Let 
us first introduce the basic ideas of the QM-cluster approach in a schematic way, 
the details and, in particular, the appropriate relativistic formulation will be given 
later. Consider a snapshot at some fixed time r. This time shold be large enough 
(> 1 fm/c), so that the "initial stage" interactions occuring in the nucleus-nucleus 
overlap zone are finished. At the given time r, we consider the locations (in R^) of 
all particles produced in the initial stage. There are, by chance, regions with high 
density and such with low density. To be quantitative, we introduce a "critical dens- 
ity" /9o, and look for domains with p > po. The high density domains, the connected 
regions with p > po^ are referred to as quark matter clusters (or droplets). Once 
they have been formed, these clusters are treated macroscopically, characterized by 
a distribution function for energy density ej(x, r), momentum pj(x, r), and fiavour 
/j(x, r), with i referring to cluster i. There are no constituents explicitly kept track 
of, and there is no memory referring to the production process. 

In order to determine high density regions at given time r, we proceed as follows: 
we assign a "critical volume" Vq = l/po to all particles. Domains of high density 
correspond to overlapping volumes (see fig. 7). So the task here is to find connected 
objects with at least pairwise overlap of individual volumes (in fig. 7, we find 
four such objects). These objects are now considered as quark matter clusters (or 
droplets). It should be noted that the "critical volume" Vq is not the usual volume 
47ri?)[/3 of a hadron, it is rather defined by the requirement that if hadronic matter 
is compressed beyond Vq per hadron, the individual hadrons cease to exist and quark 
matter is formed. So Vq is less than 47ri?|/3. Rather than the critical volume Vq or 
the critical density po, we usually use the corresponding critical energy density Cq. 
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Figure 7: Overlapping volumes, rep- 
resenting high energy densities. 



Figure 8: The collision zone (in the t- 
z plane) of an A A collision. Each dot 
represents an origin of a string evolu- 
tion. 



This critical energy density eo (or equivalently Vq or po) is the crucial parameter of 
our approach. 

For the first stage, we use the independent string model, to be more precise 
the basic VENUS model without final state interactions [3]. In fig. 8, we show a 
typical event: we consider a projection to the t-z plane, nucleon trajectories are 
represented by straight lines, interactions are represented by dots. Some nucleons 
interact (participants), some survive the interaction zone (spectators). Each dot, 
representing interaction, is a point of string formation, or, in other words, the origin 
of a string evolution. What happens after a string formation point? We use the 
standard procedure of classical relativistic string dynamics and decays [3]. In fig. 
9 a typical example of the space-time evolution is shown. The upper rectangles 
represent produced particles (hadrons and resonances), the arrows indicate particle 
trajectories. Remarkable is the strict ordering of the directions, being a consequence 
of the covariant string breaking mechanism. 

Being able to construct in a first stage, event by event, particle trajectories 
defined by their origins in space and the four-momenta, we can proceed to stage two, 
the analysis of energy densities and cluster formation at fixed time. The question 
is, what we mean by fixed time, which frame we are using. In this context, it is 
important to realize the empirically found correlation between the average rapidity 
y and the space-time rapidity C, 

y = C, (17) 

with deviations only around z k, t due to the finite size of the nuclei. So there 
is an "inner region" with Bjorken-type behaviour y = C? where all the particle 
momenta point back to the origin, and an "outer region" with parallel velocity 
vectors. Correspondingly, the dotted line in fig. 10, a hyperbola in the inner region 
and the tangents at z^ and z'^ ^ represent equal proper time r (on the average). In 
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Figure 9: Trajectories of string frag 
ments. 




->- z 



Figure 10: Space-time evolution of 
nucleus-nucleus scattering. The dot- 
ted line represents constant proper time 
(big dots: hyperbola, small dots: tan- 
gent). 



this way, obviously, also an average comoving frame is defined. The hypersurface 
defined by the hyperbola/tangents (dotted line) together with arbitrary x ^ is called 
r-hypersurface. We are now in a position to specify the frame for interactions: we 
investigate densities (or overlap) at constant r, which means on r-hypersurfaces. 

Having specified the frame and correspondingly the time coordinate r, we have 
to introduce a useful longitudinal coordinate. We use the "proper length" 



s :- 



dz* , 



:i8i 



with an integration at constant r, and with dz* being a longitudinal length in the 
average comoving frame defined by the r-hypersurface. 

To specify the geometrical properties of particles or clusters, we use the variables 
r, s, and f± = (rx^ry). At given r, a particle or cluster is considered as cylinder in 
s, r^;, Vy space, with the axis along the s-axis. The object is characterized by a lower 
and upper value of s, Si and S2, and a transverse radius r±. 

We are now in a shape to construct high density domains, to be considered 
as quark matter clusters. As discussed above, these domains are constructed, as 
in percolation models, by investigating geometrical overlaps of individual objects. 
There are two types of objects: initially, we have only hadrons and resonances, later 
also clusters contribute, which have been formed earlier. In any case, at a given time 
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r, the objects are considered cylindric in (s, r^;, rj^)-space. Connected overlapping 
objects in this space define clusters. Such a cluster has in general a very irregular 
shape, which is "smoothened" in the sense that this complicated shape is replaced 
by a cylinder of the same volume, the same length S2 — Si, and the same center 

Starting at some initial time Tq (presently 1 fm/c), we step through time as 
Tj+i = Ti + 5t^ constructing clusters at each time r^. For the time evolution of 
clusters, we presently assume purely longitudinal expansion, 

AC(r,+i) = AC(rO. (19) 

The last and most difficult topic to be discussed is the hadronization. The power 
of our percolation approach is that any hadronization scheme can be plugged into 
our approach and tested in a very detailed fashion. This is what we plan for the 
future. Currently, we present a very simple scenario, which can be implemented quite 
easily. Since our cluster expands, the energy density decreases, and drops at some 
stage below Cq. Per definition, we hadronize the cluster at this point instantaneously. 
Since the clusters turn out to have essentially longitudinal shape (As ^ Ar^), we 
proceed as follows. The cluster is cut into many short pieces in s, all of them having 
the same mass m with the requirement of m being around some parameter niseg- 
The small clusters then decay isotropically according to phase space [18]. So this is 
essentially the decay of many fireballs at different rapidities. 

Crucial for the approach is the initialization, i.e., the volume Vq assigned to the 
hadrons and resonances. The "critical" volume Vq, or equivalently the "critical" 
energy density eo, is therefore a most important parameter. One believes nowadays, 
that the QCD phase transition is not first order, but rather a narrow transition 
within a range of about 10-20 MeV, corresponding to a e-range of maybe 0.1 to 
2 GeV/fm"^. So the critical energy should be somewhere in this range, perhaps 1 
GeV/fm"^. We calculate cluster-volume distributions for values of e (CED) of 0.15, 
0.50, and 1.0 GeV/fm"^, as shown in Fig. 11. The numbers are not normalized, we 
show the number of clusters per volume bin AV, found in 1200 simulations; the 
bin sizes for the three distributions are, from top to bottom: Ay=10 fm"^, 20 fm"^, 
60 fm"^. From the grey-scales, we can also read off how the clusters are distributed 
in energy density. We show only results for r = 2 fm/c, since volumes and energy 
densities turn out to scale in a trivial manner, as V ~ r and e ~ r~^, so the 
distributions for different r look similar up to a scale transformation. 

We observe the following: for eo = 1 GeV/fm"^, the distribution peaks at small 
values of V, dropping very fast with increasing volume V. Reducing the CFD to 0.50 
GeV/fm"^, the distribution gets wider, roughly by a factor of two. Although, as for 
£0 = 1 GeV/fm"^, very small clusters are favoured, the fiuctuations are considerably 
larger. Reducing eo further to 0.15 GeV/fm"^ (= nuclear matter density), we observe 
a drastic change in the distribution, a maximum at large values of V emerged. This 
means, it is quite likely to observe just one big cluster (filling the whole volume). 
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Figure 11: Distribution 
of cluster volumes for dif- 
ferent values of the critical 
energy density (CED) to- 



Considering Eq = 1 GeV/fm"^ as the most realistic value, the upper plot of fig. 
11 represents the realistic world (the lower plots are just mathematical exercise). 
Although here it is most likely to produce just small clusters (mainly hadrons and 
resonances), there is nevertheless a reasonable probability to form big clusters, with 
our statistics of 1200 simulations up to 65 fm"^. The question is, whether such 
"miniplasmas" can be "isolated" in event-by-event experiments. 



5 Conclusions 



We discussed the basic concepts of the string model, also to demonstrate that it is 
theoretically well justified. We outlined the limitations of the independent string 
model at low energies, at high energies, and at high particle densities. In this context, 
we reported on recent attempts to overcome these limitations by combining soft and 
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hard scattering and by considering secondary interactions. 
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